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A simplified model of quantum electrodynamics involving a charged two-state sys- 
tem interacting with an electromagnetic field mode is examined. By extending the 
Schrodinger equation to include stochastic and nonlinear terms the dynamical pro- 
' cess of quantum state reduction can be represented. A specific choice of modified 
' Schrodinger dynamics is shown to result in stable coherent field states. The two-state 
system undergoes an induced state reduction to a generalised current state due to 
, its interaction with the field mode. Numerical results are presented demonstrating 
I state reduction dynamics for an initial superposition of two current states. An in- 
^ I duced reduction time-scale for the two-state system is derived and confirmed by the 
' numerics. 

Q^. PACS numbers: 03.65.Ta, 12.20.-m 

O^. I. INTRODUCTION 
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Standard quantum dynamics cannot describe how a particular measured outcome is re- 
alised in an experiment — this feature underlies our failure to unify classical physics with its 
G^ . quantum roots. A direct response to this problem is to take the view that standard quan- 
^ ! tum theory is incomplete and that a complete quantum theory should provide a dynamical 
^ I description of quantum state reduction. In this way it is hoped to provide an objective 
Q ■ description of quantum measurement without recourse to a classical measuring device. 
OO , The great success of quantum theory does not leave much room to manoeuvre when 

^ ' proposing alternatives. Any new theory must agree with the experimentally verified pre- 
dictions of quantum theory. This roughly implies that the dynamics should reduce to 
^ ■ Schrodinger dynamics for micro systems and that state reduction effects should become 
more important with system size. More specifically, the form of any modifications we make 
should involve a random factor (at some level of approximation) to reflect the random na- 
ture of state reduction, and we would expect some nonlinearity in the dynamics, reflecting 
a feedback from the state vector to the probability of the outcome [if. 

These features can be adequately represented by a stochastic differential equation for 
the quantum state [2I, 0] (for general reviews see 0,111), the most well known model of this 
type being the continuous spontaneous localisation (CSL) model [oj. In the CSL model, 
state reduction occurs onto a basis of smeared position states leading to definite localised 
positions for systems of macroscopic size. However, there is scope for further model build- 
ing, principally because we can choose from many other preferred bases for state reduction 
(e.g. 0, [Si, i9i])- Until we have observational evidence of state reduction dynamics we are free 
to experiment with different proposals. 
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Here we investigate the state reduction process for a simplified form of quantum electro- 
dynamics. We consider a charged two-state quantum system interacting with a quantum 
electromagnetic field mode. We investigate a model which describes state reduction onto a 
coherent field state basis. State reduction for the two-state system happens indirectly via 
standard interaction with the field mode. We find that the two-state system evolves towards 
a generalised current eigenstate. 

The appealing feature of this model is that the stable eigenstates of the system are a 
natural description of the classical world whereby the field is represented by a coherent 
state and the two-state system (representative of charged matter) by a current eigenstate. 
Modifications to Schrodinger's equation only involve the field mode — it is not necessary to 
make further modifications in order to independently describe state reduction in the two- 
state system. Also, it is possible to generalise the model to field theory without generating 
infinite particle creation [l^ (this is a problem which appears, e.g., with the CSL model 



111 . |12|). although it should be pointed out that this results in a loss of Lorentz invariance. 



The idea that dynamical state reduction might be induced from one field to another is 
not new. It has appeared in previous attempts to describe state reduction in fermionic field 
theory by invoking a scalar field with state reduction properties to which the fermion is 
coupled 11, 13]. To our knowledge the process has not been previously studied in numerical 



detail. 

In the next section we outline our model, detailing the Hamiltonian for our simplified QED 
system and the modified stochastic dynamics of the state vector. In section IIIII we consider 
the evolution of a specific initial state vector composed of a superposition of current states. 
We outline the state reduction process and derive a characteristic time-scale. In section [IV] 
we present some numerical results which confirm the analytic results and intuition gained 
in section IIIII We end with some conclusions. 



II. THE MODEL 
A. QED Hamiltonian 

Our simplified QED Hamiltonian is given by 

H = ua^a + \i'crz + go'xia + a^). (1) 

The (TjS are the standard Pauli matrices and a, are the annihilation and creation operators 
for a standard harmonic oscillator. The parameters u, u, and g are constants which we 
interpret as the field mode frequency, the two-state energy gap, and the coupling respectively. 
Models of this type are commonly used in quantum optics to describe the interaction of 



a two-state atom with an electromagnetic field [13|. In particular, a further approximation 
to the interaction term to exclude processes such as those where a photon is created as the 
atom is exited (the rotating wave approximation), leads to the celebrated Jaynes-Cummings 
model. Alternatively, this model can be derived from full QED by restricting space to just 
two points (time remains continuous) and considering only one photon polarisation. The 
result is two fermion states and one field mode. This simplification is inadequate for a 
number of reasons. Most importantly, the model can no longer be described as a gauge 
theory as there is no local symmetry associated with the fermion field. We also lose the 
concepts of conserved 4-current and relativistic invariance. 
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On the other hand, our simphfied Hamiltonian has the structure of the QED Hamiltonian: 
it is composed of free matter and field parts and an interaction between a current operator 

and the field. (That corresponds to the current can be seen as we restrict full QED 
to a space of two points.) These features suffice in order to make a tentative attempt to 
understand quantum state reduction in QED. 

Since this model can be applied to quantum optical systems, it is worth pointing out that 
stochastic adjustments to Schrodinger dynamics have been used in the past to describe the 
effect of an environment on optical and atomic behaviour 
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B. Stochastic Schrodinger equation 

We shall consider a modified Schrodinger equation of the form 

dlV^t) = I [--iH - X^a^a + }?a\a - \\^a\a^ dt + A (a - dBt - iAaJdE^j (2) 

where A is a constant coupling parameter with units of [time]~^/^ and 

at = {^pt\a\iJt)- (3) 



The complex Wiener process Bt satisfies dBtdB^ = 2dt, and dBtdBt = dB^dB^ = [Ij 
Note that the Schrodinger equation is recovered by setting A = and that the additional 
A-dependent terms only involve field operators. 

Equation ([2]) is a stochastic differential equation for the state. It can be derived by looking 
for a state diffusion process describing stochastic fluctuations in the quantum probability 
amplitudes of possible a-eigenstates [4;]. The probabilities for the stochastic fluctuations must 
be given by the quantum probabilities at each point in time. This ensures that stochastic 
outcomes match with the initial quantum predictions (and leads to nonlinearity) . The 
equation has the property that general states evolve towards a-eigenstates on a time-scale 
approximately inversely proportional to the variance in a (see below) . 

Using the rules of Ito calculus we can show that 

d{^tm = o, (4) 

i.e. the norm of the state is preserved under the dynamics ([2]). This allows us to maintain 
the quantum interpretation of the state as a probability amphtude. 
The process for the conditional expectation of a is found to be 

da* = {-luat - ig{a,)t - X^at] dt + X{{Aaf)tdB + X{\Aa\^)tdB\ (5) 

where Aa = a — at and we have used the convenient notation {■)t = {ipt \ ■ \'4't)- When the 
state \ipt) is an a-eigenstate, the process at is nonstochastic. In this case, the effect of the 
modiflcations to Schrodinger's equation is a decay in the magnitude of at at a rate given by 
A^. This can be interpreted as being due to the fact that the classical stochastic process Bt 
in equation ([2]) is only coupled to photon annihilation events. This results in a gradual loss 
of energy. Appropriate tuning of A is required to ensure that fleld energy loss is negligible 
whilst state reduction timescales are appropriately small for strong flelds 13 . 



In order to demonstrate the state reducing properties of equation ([2]) we consider the 
conditional variance in a deflned by 

Vari(a) = {\Aa\'')t (6) 
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The conditional variance process is found to be 

dVart(a) = ig {{a^{a - a^))t - {cr^)t{ia - a''))t} dt 

-2X^Vait{a)dt - {|((Aa)')t|2 + Var2(a)} dt 
+X{\Aa\'^Aa)tdBt + X{Aa^Aa\'^)tdBl (7) 

The first hne on the right side is the effect on Varf(a) due to the interaction between the 
two-state system and field. We see that Vart(a) increases if there is a negative conditional 
covariance between ax and —i{a — a^). As Vart(a) increases, other terms on the right side 
of equation ([7]) become more important. 

If we ignore the interaction with the two-state system the variance process is a super- 
martingale. Integrating and taking the unconditional expectation of equation ([7]) we then 
find 

E[Var((a)] = Varo(a) - 2A^ [ dME[Var„(a)] - 2A^ [ duE [\{{Aaf)u\'^ + Var^(a)] . (8) 

Jo Jo 

This demonstrates that the expected variance will decrease in time and therefore that the 
field approaches an a-eigenstate. 

We can estimate the time-scale Tq for state reduction to an a-eigenstate from the following: 

E[Vart(a)] - Varo(a) t 
Varo(a) 

By freezing the stochastic terms on the right side of equation ([8]) at time t = we find 

Varo(a) 

'"'^^ A2{|((Aa)2)o|2 + Varg(a)} 

Even though A must be very small to prevent energy loss, the reduction time-scale can 
be small if the variance in a is large. This demonstrates the specific way in which the 
effectiveness of state reduction is dependent on the size of the system. A superposition of 
greatly differing a-eigenstates has large variance and state reduction will occur rapidly; if 
Varo(a) ~ |((Aa)^)o| ~ 0{1) then state reduction takes an extremely long time due to the 
small value of A. In principle, the dependence of the reduction time on the a-variance is a 
testable result. An experiment designed to test this might look for diminishing quantum 
interference as a signature of the state reduction [T^. The technical difficulty would be to 
eliminate any environmental influence. 

In order to see that stochastic probabilities for outcomes match with the initial quantum 
estimates we consider an initial superposition of coherent states 

l^o) = X^Ciltti), (11) 

i 

where it is supposed that the coherent states \ai), where a\ai) = ai\ai), are sufficiently 
separated in phase space that {ai\aj) ~ 6ij (this requires that \ai — a^l 1 for i 7^ j). We 
assume for simplicity that H = since we wish to focus on the reduction dynamics. 

Note that the coherent states are eigenstates of all terms in the evolution equation ([2]) 
except the term —X^a^a. This means that most terms in the evolution equation simply mod- 
ify the coefficients q, leaving the set of states {jai)} intact. The term — A^a^a is responsible 
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(10) 
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for diffusing the initial set of coherent states. In order to analyse its effects, consider a state 
evolution involving just this term: 

^IV.,) = -W^a\^t). (12) 

This equation has the solution 

= 5^c,exp|i(|a,e-^'f - \a,\')] la.e"^'*). (13) 

i 

This indicates that the individual coherent state contributions to the overall superposition 
will decay to the ground state on a timescale A~^. For times t -C we can ignore this 
effect and assume that the set of coherent states contributing to the superposition is fixed 
as From equation (fTOl) we see that ^ for large a-variance. We can therefore 

assume that in the region Ta <^t <^ \~^, the state has reduced to one of the coherent states 

We can define an approximate coherent state projection operator by = \ak){<yk\- 
For Ta t <^ A^^ the conditional expectation (iptlPaJipt) is approximately either 1 or 
depending on whether the state is reduced to \ak) or not. From equation ([2]) we find 

dii^tlP^J^Pt)} = Mokit)\\ak - at)dBt + X\ckit)\\al - aJ)dS;, (14) 
from which we conclude that 

(HPa.m = E[(^i|P„J^<)|^o] = E[l|v,) = |a,)|^0], (15) 

valid for ^ t ^ A^^. The left hand equality gives the initial quantum probability for the 
outcome \ipt) = the right hand equality gives the stochastic probability for the same 
outcome. 



III. AN EXPERIMENT 

Consider an initial situation where the field is in a coherent state |a) and the two-state 
system is in a superposition of current (ax) eigenstates: 

l^(0)) = {ci(j)+C2(^_\)}|a). (16) 

In order to understand the behaviour of this model without explicitly solving it we consider 
some possible scenarios depending on the relative sizes of the model parameters. 

Since the field state is initially coherent, the variance in a is initially small. As this is our 
measure of system size, we can assume that equation ([2]) initially reduces to the Schrodinger 
equation. 

A. CO, V domination 

First consider the case where the coupling g tends to zero. We then have 
\i){t))g=Q = exp {-it (cua^a + iz/a^)} IV^(O)) 

= {ci(t)(j)+C2(t)(_\)}|e-*-«), (17) 
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where 

Ci(t) = Ci cos |tz/ — ZC2 sin |tz/, 

C2(t) = C2 cos |tz/ — ici sin itz/. (18) 
Here we find independent state-space rotations for the two-state system and the field. 



B. g domination 



Next consider the case where the field frequency uj and the two-state energy gap u tend 
to zero. In this limit we find the following solution to the Schrodinger equation 

\ip{t))^=^=o = exp{-itgax{a + a')}\ip{0)) 

= c,e'^^ ( I ) l(« - ^9t)) + c^e'^' ( _\ ) l(« + Wt)). (19) 

where 0j are real-valued phase factors. The interaction between two-state system and field 
leads to a correlation between current eigenstates and field coherent states. 

These two limits provide an overall picture of two competing processes. When the inter- 
action g is dominant, different current eigenstates become associated with different coherent 
field states. When the interaction is negligible we find that current states rotate into one 
another. 



C. A domination 

Consider the scenario uj ^ v <^ g. From equation (1191) we can estimate the growth of 
variance in a: 

Var,(a)=. Vt2|cinc2p^|((Aa)2),|, (20) 

and the conditional covariance between ax and —i[a — a^): 

CoYt{ax,-i{a- a))) = -i{ax{a - a^))t + i{ax)t{{a - a)))t 

~ -%gt\ciW2? (21) 

(the approximation results from the fact that we have assumed that different coherent states 
are nearly orthogonal). As the variance grows, A-effects come to dominate the evolution. 
This can be seen from equations (1201) and fl^Tl) by noting that the Var^(a) and |((Aa)^)tp 
terms in equation ([7]) will soon outgrow the Cov^ [ax, —i{a — a))) term. Let us define an 
approximate time s at which a transition from g to A-dominated behaviour happens. Once 
in the A-dominated phase, the a-variance will tend to zero on a time-scale given by equation 
(fTOj) . The stochastic evolution effectively chooses between the possible states (1, iy\a — igt) 
and (1, — l)^|a -|- igt) leading towards a definite field coherent state and a definite current 
state. 

Using (ITOl) and (!20!) we find that a total state reduction time-scale for the current state 
is given by 
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FIG. 1: Sample path showing growth and subsequent decay of Vart(a). The inital state involves an 
equal superposition of two different current states. The coupling between current and field mode leads 
at first to a growth in Yait{a) via standard Hamiltonian dynamics. When Var((a) is sufficiently large, 
the stochastic dynamics cause a decay in variance as the a-state collapses. Parameters are uj = u = 0.5, 
g = 4.0, A = 0.2. 



In order to determine s we find the value which minimises To-^ . Setting 

d 1 

= ^> 1 ~ 

ds A^(7^s^|cip|c2 

therefore 



~ p/3^2/3|c^|2/3|c2|2/3- (^4) 

This is the characteristic time-scale for state reduction of a superposition of current eigen- 
states in a two-state system coupled to a field mode. 



IV. NUMERICAL RESULTS 

In order to test these results we have implemented a numerical simulator. Starting from 
an initial state of the form given in equation (fT6|) with coefficients Ci = C2 = 1/2 and with 
a = 4, we have evolved the state according to equation ([2]). Figures 1 and 2 show sample 
paths for Vart(a) and {crx)t respectively. Here we have chosen = z/ = 0.5, g = 4.0, A = 0.2. 
This choice of parameters are such that g dominates the initial evolution. The parameter 
A is chosen to be large enough to observe state reduction for the modestly sized numerical 
system, but small enough that energy loss is not too significant. 

It is seen that the a-variance initially grows smoothly, approximately proportional to 
(see equation fl20|) ). After a certain time the evolution becomes dominated by stochastic 
movements with a notable downward drift. The a-variance finally becomes stable again as 
its value tends to zero (at around t = 1.5). 
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FIG. 2: Evolution of the expected current {ax)t- An initial (equal) superposition of two different current 
states {{(Jx)o = 0) undergoes state reduction to a definite current state {\{(Tx)t\ = !)■ Same random 
sample as figure 1. 



The current in figure 2 corresponds to the same sample. It is seen that the two-state 
system tends stochastically towards a definite current state — in this particular case (1, 1)^ 
This happens at the same time as the a-variance goes to zero. 




FIG. 3: Expected current reduction time as a function of the current-field coupling g. Other parameters 
are a; = = 0, A = 0.2. The filled line shows a fit to the curve g~^^^- The error bars show the one 
standard deviation region for the reduction time estimated from the simulated data. 



We define the reduction time in our numerics as the stopping time at which | {crx)t\ > 0.99, 
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tau_{sigma_x} 

FIG. 4: Histogram of |(cTa;)t| > 0.99 stopping times for g = 5 {uj = u = 0, X = 0.2). The smoothed 
line is the Gaussian kernel density estimate. 



I.e. 

r,, = inf{t >0;|K)t| >0.99}. (25) 

Figure 3 shows numerical estimates of the expectation of this stopping time for a range 
of g values. Other parameters are chosen to be A = 0.2 and u) = u = 0. Expectations 
are determined using 100 sample paths for each point. The estimates are seen to almost 
perfectly fit the curve kg"'^/'^ with k = 3.593. This confirms the form of the characteristic 
state reduction time given in equation flMl) . The error bars indicate a range of one standard 
deviation for . We observe that the variance in the stopping time decreases as g increases. 

Figure 4 shows a histogram of t^^ results for 100 sample paths with g = 3. Also shown 
is the Gaussian kernel density estimate 3]. The distribution of stopping times is approx- 



imately normal with a slight skew. The important result is that the density of events with 
large deviation from the mean stopping time falls away rapidly. 

We find that occurrences of {ax)t — ^ 1 and {ax)t — ^ —1 are even to within statistical 
error. This agrees with standard quantum predictions for a measurement of ax given the 
initial conditions. 



V. CONCLUSIONS 

We have investigated the state reduction dynamics described by a modified Schrodinger 
equation for a simplified QED system. Our simplified QED system is composed of a charged 
two-state system interacting with an electromagnetic field mode. The modified dynamics 
result in state reduction to an over-complete coherent field state basis. At the same time, the 
electromagnetic interaction leads to the development of correlation between current states 
and field coherent states. This leads to an induced state reduction to a current state basis 
for the two-state system. We have considered evolution from a specific initial condition 
composed of an equal superpostion of the two possible current eigenstates with a coherent 
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field state. We have implemented a numerical simulator to generate sample state evolutions. 
These have confirmed our analytic predictions for the general behaviour and the time-scales 
involved. 

The time-scale for state reduction of the the two-state system is found to behave as 
g-'^/^Yo^-'^l^ (cTx) . We might speculatively suggest that the reduction time-scale for more 
general charged matter systems might scale as Var~^/^(j^) for 4-current j^. 
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